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Abstract. We propose a method to extract nnclear charge distribntions from 
elastic electron scattering data based npon a mean held approach. The nnclear charge 
distributions are generated by solving the Schrodinger equation with a mean-held potential 
expanded in terms of Hermite functions. The coefficients of the potential are changed in 
order to obtain the best ht with the cross section data. Applications to the ^^(7, 
and nuclei are presented. 
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1. Introduction 


The distribution of the electromagnetic charge in the nuclear ground state is one of the 
more interesting, and investigated, properties of the atomic nuclei. Nowadays, the widest 
and most accurate set of information about this quantity is provided by elastic electron 
scattering experiments. Ultra relativistic electrons are scattered by the electrostatic po¬ 
tential generated by the nuclear charge distribution. From the measured cross section one 
extracts information about this distribution. 

The straightforward approach to get the charge distribution is the solution of the 
inverse scattering problem. Unfortunately there are principle and pragmatical difficulties 
making this approach extremely impractical [1] and for this reason it has never been applied 
in this context. 

The approach which has been commonly used consists in solving directly the scattering 
problem making a guess about the form of the potential. The parameters hxing the form 
of the potential are changed until the experimental cross sections are reproduced. 

In early times, the quality of the data available was such that they could be htted 
reasonably well using simple functions containing few parameters (e.g. Fermi distribu¬ 
tions). In the last two decades, thanks to the the development of experimental facilities, 
there has been a big improvement of both quality and quantity of the data. Reasonable 
hts of modern elastic cross section data can be obtained only increasing the number of 
the free parameters. Instead of inserting a dependence of new parameters in the old func¬ 
tions describing the densities it results more convenient to make use of the so-called model 
independent techniques [2,3]. 

These techniques consist in expanding the charge distribution on a orthonormal basis 
and changing the coefficients of the expansion to achieve the best ht of the cross section. 
In principle the model independent techniques are very general and clean, but their ap¬ 
plication to the specihc problem of htting experimental cross sections shows ambiguities 
and uncertainties. The expansion of the density on a basis does not converge but rather, 
with the increase of the number of the expansion terms shows an enhancement of the 
uncertainty band of the charge distribution. All this happens in spite of the fact that the 
quality of the ht to the cross section remains constant. We have shown in ref. [4] that this 
instability of the method is not related to a specihc choice of the expansion basis, but it is 
rather a general drawback of the model independent techniques. 

Pragmatical recipes have been developed to keep under control the instabilities of 
the model independent techniques. These are ranging from analytical continuations of the 
experimental data in the unmeasured region [3] to the identihcation of the optimal number 
of expansion coefficients [4]. 

This unsatisfactory situation still rehects the fact that it is not possible to solve the 
inverse scattering problem for a set of data known on a discrete and hnite grid and with 
statistical error [1]. 

In the present article we propose an extraction method alternative to the model inde¬ 
pendent techniques. This method is in principle less general, but it produces more stable 
solutions. 

Our idea is to use the mean held model of the nucleus to extract the charge distribution 
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from the experimental data. The nucleus is described in the framework of a conventional 
shell model. The charge distribution is given by the convolution of the nucleon electro¬ 
magnetic form factor with nucleon densities generated by the sum of the squares of the 
single particle wave functions. The parameters describing the nuclear average potential 
generating the wave functions can be changed to ht the cross section data. 

It is clear that this approach has more physical insight than the model independent 
techniques, and it is restricting the class of functions used for the description of the charge 
distributions. For example in this case the charge distribution is ensured to be always 
positive; this is not the case for the model independent techniques. 

For different purposes, hts of elastic electron scattering data using the nuclear shell 
model with an energy dependent mean-held, have been published already in the 1967 [5]. 
The basic idea of this work was to use the experimental data in order to hx an empirically 
well grounded mean-held potential. 

For us the shell model is only a tool which imposes reasonable physics constraints 
on the class of functions used to build up the charge distributions. We are showing that 
energy independent mean-helds with enough degrees of freedom can reproduce elastic 
electron scattering data on a much wider momentum transfer range than that available in 
1967. In spite of this we are not attributing any special meaning to the potentials obtained. 

In addition, we do not consider the charge distributions we present in this work as 
the most realistic ones. We are conscious that a careful extraction of the nuclear charge 
distribution should take care that the scattering data are affected by dispersion effects, 
meson exchange currents, neutron charge distribution and other phenomena which are 
ignored in this work. 

The aim of the present paper is to present an extraction method providing charge 
distributions which result to be much more stable against the number of the free coefficients 
than the traditional model independent techniques. 

In the next section we describe the method used to generate the mean held potential 
and we test its performances on pseudodata. The method is applied in sect. 3 to reproduce 
the experimental cross section of the ^^(7, ^^0, and nuclei. 


2. A generalized mean-field model. 

As discussed in the introduction, we are working in the framework of the direct scat¬ 
tering approach. We build up a charge distribution Pcir) generating the electrostatic 
potential used to solve the Dirac equation for ultra-relativistic electrons. The techniques 
used to solve the Dirac equation, and to obtain the elastic scattering cross section, are 
presented in ref. [6]. 

The charge distribution is obtained folding the nucleon density distribution p(r) with 
the nucleon form factor /n(r',r): 


Pc(r) 


= j pW)fnW,v)d^v' 


( 2 . 1 ) 
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The density distribution of a system of A particles is defined as: 


A 

p(r) =< Tl ^5(r — ri)|\l/>, (2.2) 

i=l 

where we have indicated with |\1/ > the wave function describing the ground state of the 
system. 

In the mean field approximation |\1/ > is a Slater determinant of single particle wave 
functions (/>a(r) and the density can be written as: 


p(^) =^<P»(r)*<Pa(r), 

a 


(2.3) 


where a indicates the set of quantum numbers identifying the single particle level. 
We work with a mean field potential of the form: 

VMF(r) = V(r) + -^V(r) 1 ■ s + Vcoui(r) (2.4) 

r dr 

where Vcoui indicates the Coulomb potential which has been chosen as the electrostatic 
potential generated by a uniformly charge distribution of radius 1.2^3 fm. 

In this expression we should fix the expression of the central term V (r) and the 
strength A of the spin-orbit term. 

Our method consists in expanding V{r) on a basis of orthonormal functions Pni^)'- 


OO 

V{r) = '^CnPn{r). (2.5) 

n=0 

Specifically we have chosen Pn{r) = Un{r/(3) where the Un are Hermite functions 
defined as [6]: 

Un{x) = ( 2 . 6 ) 

where we have indicated with the Hermite polynomial of order n. The orthogonality 

relation of the Un functions is given by: 

/ OO 

Un{x)Um{x) dx = 5n,m- (2.7) 

-OO 

The calculations we have performed consist in solving the single-particle Schrodinger 
equation with the potential given by the eqs. (2.4) and (2.5). With the single particle wave 
functions obtained in this way we construct the density distribution, eq. (2.3), and then 
the charge distribution, eq. (2.1). With this charge distribution we build up the Coulomb 
potential to solve the Dirac equation and to obtain the cross section which is compared 
with the experimental data. The procedure is repeated modifying the parameters Cn of 
the expansion of the potential until the minimum value of the is reached. 
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As a general rule we have populated the single particle levels which are usually thought 
to be occupied in the conventional shell model description of nuclear ground states. 

We have tested the procedure described above applying the method to reproduce 
^‘^^Pb pseudodata generated by a mean-held model with Woods-Saxon potential. The 
pseudodata have been obtained substituting in eq. (2.4): 


V{r) 


-Vo 

1 -t- exp{'^-^) 


( 2 . 8 ) 


with the parameters given in Tab.l. 

We have tested the sensitivity of our results to the choice of the nucleon form factor, 
comparing charge distributions and cross sections calculated with the same single particle 
wave functions but with various form factors. The differences between the various results 
are really negligible, contrary to what happens in the quasi-elastic peak region [8]. We 
have used the form factors given in the reference [9], but neglecting the small contribution 
of the neutrons. This means that the sum in equation (2.3) runs only on the protons. 

The charge distribution generated by the Woods-Saxon potential has been used to 
produce cross section pseudodata under the same kinematical conditions of the 
experimental data of ref. [10,11]. The percentage errors on the cross section pseudodata 
have been taken equal to the experimental ones [10,11]. 

We found that the straightforward application of our minimization techniques, based 
on the gradient method, had convergence problems. This was due to the fact during 
the variation of the parameters, some of the populated single particle levels go in the 
continuum, producing discontinuities in the charge distribution and therefore in the x^. 

In order to force the parameters to move only in regions in which all the conventional 
shell model single particle states are bound, we applied the gradient minimization technique 
to the sum of the x^ of the cross section and of the x^ of the single particle levels. 

To the above sum we added a penalty term: 


P = (2.9) 

i 

where the index i runs over all the points of the grid used in the numerical integration of 
the Schrodinger equation where V{ri) > 0. 

We have added this term to insert a reasonable physics input in order to speed up 
the minimisation procedure guiding the algorithm towards regions, in the space of the 
parameters, physically meaningful. This is the analogue in a model independent techniques 
of including a penalty for negative charge distributions. In all the hts performed, the hnal 
value of the penalty term resulted to be compatible with zero. 

In this way the merit function which is minimized depends on the strength of the 
penalty parameter p and on the errors entering in the x^ associated to the single particles 
levels. Simple recipes hxing these quantities do not exist. Our experience indicates that 
only repeated minimization performed by varying the relative importance of the terms of 
our merit function produce satisfactorily results. 

Using these prescriptions we obtained excellent hts to both the cross section and the 
single particle energies. 
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Some results of the fit to the pseudodata are shown in fig. 1. In the upper panel we 
compare the density obtained by the fit to the cross section with 8 Hermite parameters, 
with the original Woods-Saxon density shown by the full line. On the scale of the figure 
the lines are overlapped. In the other two panels of the figure we compare the central and 
the spin-orbit part of the Woods-Saxon potential (full lines) with the potential generated 
by the Hermite expansion given in eq (2.4) (dashed lines). 

The central part of the potential is rather well reproduced, especially on the surface 
region. In the internal region the potential generated by the Hermite expansion shows 
small oscillations. 

While the central potentials are rather similar, the spin-orbit potentials show big 
differences. The spin-orbit potential generated by the Hermite expansion (dashed line) 
has big oscillations in the interior of the nucleus (for values of the radius smaller than 
5 fm). This behavior is not present in the Woods-Saxon spin-orbit term, which, in this 
region, is smoothly going to zero. On the other hand, the two potentials are rather similar 
in the surface region. This behavior is not surprising but it comes out from the fact that 
we chose to determine the spin-orbit potential from the derivative of the central term [12]. 
As we have discussed before the central term are similar on the surface but they differ in 
the interior by the presence of oscillations which are generating the curious behavior of the 
spin-orbit term in this region. 

To get rid of the oscillations in the internal region of the the spin-orbit term, we have 
repeated the fit constructing this term in a different ma nn er. We calculate the derivative 
of the central potential starting from the external part of the nucleus, and after the first 
maximum of the derivative we build up the rest of the spin-orbit potential as symmetrical 
f un ction with respect to this point. The dashed-dotted lines of fig. 2 have been generated 
with this method. 

The results obtained in this way are very similar to the previous ones. This tells 
that the internal part of the spin-orbit potential has very little influence on the results. 
This fact can be understood considering the fact that only the s-waves are sensitive to the 
details of the potential in the interior of the nucleus. On the other hand the spin-orbit 
potential does not act on these waves. 

In the model independent procedures the increase of the number of expansion terms 
is producing a widening of the uncertainty band of the charge distribution. We have 
tested the stability of our model against the increase of the number of the expansion terms 
repeating the fit to the cross section with different number of expansion terms. The results 
of this study are shown in fig. 2. One can see that there is a remarkable stability of charge 
distribution. This stability of the charge distribution does not correspond to a stability of 
the potential, which for a large number of coefficients starts to have an oscillating behavior. 

These oscillations are produced by the high components of the expansion (2.5). This 
effect is rather well known in the model independent extraction of the density [2,3,4] and 
it is the origin of all the problems connected with the increasing of the uncertainty of 
the charge distribution. What fig. 2 is showing it is one of the advantages of performing 
a fit of the cross section using a mean field approach instead than a direct expansion of 
the density. The high components of the expansion are acting on the potential, but the 
density, generated by the sum of the single particle wave functions, remains stable. 
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3. Application to specific cases. 


The results presented in the previous section show that the mean held method we have 
proposed is able to reproduce the set of pseudodata generated by a mean held potential 
with great accuracy. We do not expect the same kind of performances when the method 
is applied to ht the experimental data. 

We should remark that the model independent expansion techniques of the density 
produce hts to the pseudodata of the same quality of that obtained by our mean held 
approach. 

In this section we present the results we have obtained using the mean held method 
to reproduce the electron scattering elastic cross sections measured on the ^^(7, 
and 208 p^ nuclei [10,11,13714,15]- 

As we have discussed in the case of pseudodata, our minimization procedure needs 
the knowledge of the single particle energies. The value of the experimental single particle 
energies we have used are presented in tab. 2 together with our results. We remark that 
the use of the levels presented in tab. 2 is sufficient to stabilize the minimization procedure. 

The agreement between the calculated energies and experimental ones is reasonable, 
but not as good as in the case of the pseudodata. The experimental values of the single par¬ 
ticle energies we have used have been obtained with diherent methods (extracted from the 
binding energies of neighbor nuclei, from the centroid energies in knock-out experiments, 
etc.). We are identifying these energy values with the eigenvalues of the single particle 
hamiltonians building up the mean held many body hamiltonian. This identihcation is 
quite extreme and doubtful. Many-body effects present both in the ground state proper¬ 
ties and in the extraction mechanism can strongly modify the interpretation in terms of 
mean held approach. 

A comparison of the performances of our approach against those of the model inde¬ 
pendent techniques is shown in hg. 3 where the charge distributions obtained htting the 
data [11] with a Fourier-Bessel (FB) expansion of the density and with our method 
are shown. 

The dark areas represent the region dehned by the upper and lower envelopes of the 
densities and potentials obtained by a set of 100 hts of the experimental cross section 
normally sampled within the error band [4]. The FB densities, shown in the panels A and 
B, have been obtained with 9 and 13 coefficients respectively. In the panels C and D we 
show the densities obtained with our model and the mean held potentials (right vertical 
scales and lighter areas). In C the ht has been done with 8 expansion coefficients, while in 
D with 14 coefficients. 

The distribution bands of the densities and the averaged reduced show that the 
quality of the hts obtained with our model is analogous to that obtained by the usual 
model independent techniques in their better performances. In the FB case the increase of 
the number of free parameters generates large charge distribution bands. With our method 
the instability generated by increasing the number of the expansion coefficients shows up 
in the potential, but the charge distributions remain quite stable. 

In hg. 4 we compare the results of our method with those obtained using the FB 
expansion of the density for all the nuclei we have considered. In this hgure the full lines 
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have been obtained with the shell model, while the dispersion band has been obtained by 
a set of FB hts. As we have seen in hg. 3 the dispersion band of onr model is of the same 
order of magnitnde of that of the FB hts. 

The fnll lines are completely overlapped by the FB density bands. Only in the central 
region of we notice some difference between the FB density band and the fnll line 
prodnced by the shell model ht. 

In tab. 3 we compare the valnes of the rednced obtained for varions hts. We 
observe that the mean held hts have valnes comparable to those obtained with the FB 
expansion and slightly larger than those obtained with the Hermite expansion. 

We do not insert hgnres of the cross sections, since on the nsnal printable scale, the 
experimental data and all the calcnlated lines are overlapping, as clearly indicated by the 
valnes of the of tab.3. 

The potentials generating the single particle wavefnnctions which prodnce the densities 
htting the experimental data are shown in hg. 5. We observe the odd behavior of the spin- 
orbit terms which are generated as derivative of the central terms. The central terms have 
qnite a few oscillations and this generates big variations in the spin-orbit term. We have 
seen in the previons section that the resnlts are not very sensitive to the details of the 
spin-orbit term, however the shape of the potentials in the lower panel of hg. 5 is qnite 
odd. 

To get rid of this problem we have snbstitnted the spin-orbit term given in eq. (2.4) 
as derivative of the central term, with a Woods-Saxon term of the form (2.8). The three 
additional parameters of the spin-orbit term are also entering in the htting procednre. The 
resnlts of these hts are presented in tab. 2 and 3 where we have called these resnlts as 
SMHWS (Shell Model Hermite Woods-Saxon). 

The central potentials of the Hermite model and of the Hermite pins Woods-Saxon 
are compared in hg. 6. One can see that the SMHWS potentials are smoother than the 
Hermite potentials, bnt in general the shapes are qnite similar. 

The SMHWS potentials have reasonable behavior from the physics point of view. We 
shonld remark that they are always negative and they are smoothly going to zero for large 
valnes of the radins. The only exception is the potential of which shows a strange 
pocket aronnd 9 fm. All onr attempts to get rid of this pocket have prodnced hts of mnch 
worse qnality. 

The two mean held models prodnce hts of the same qnality, as well as it is possible to 
see from the resnlts shown in tab. 2 and 3 and from hg. 4 where the charge distribntions 
obtained with the SMHWS approach are presented by the dashed lines. 

The resnlt obtained so far show that the mean held approach prodnces hts of the same 
qnality of those obtained nsing the model independent techniqnes. In hg. 6 we present the 
densities obtained by onr hts in logarithmic scale. The densities prodnced by the model 
independent techniqnes (bot with FB or Hermite expansion basis) show clear anomalies in 
the tail, while the mean held densities have the proper exponential decay compatible with 
the qnantnm mechanics reqnirements. 


3. Conclusions. 



In this paper we have presented a method, alternative to the model independent 
techniqnes, to extract the charge density distribntions from elastic electron scattering 
data. 

The method is based on the mean held model of the nnclens. The mean held potential 
is expanded on a orthonormal basis of Hermite fnnctions, and the expansion coefficients 
are modihed in order to reprodnce the cross section. 

Clearly this method is restricting the set of fnnctions describing the charge distri- 
bntion, while model independent techniqnes are more general. On the other hand, this 
restriction is stabilizing the solntion of the problem. We have shown in sect.2 and in hg. 
3 that increasing the nnmber of coefficients creates oscillation in the mean held potential, 
bnt these oscillations are not prodncing instabilities in the charge distribntions. 

The loss in generality is slightly penalizing the qnality of the of the ht. In tab. 3 we 
observe that the valnes prodnced by a model independent ht with a basis of Hermite 
fnnctions are a slightly better than the valnes obtained with onr method. 

The better performances of the model independent techniqnes have the serions draw¬ 
back that the charge distribntions prodnced have an nnphysical behavior at high valnes 
of the nnclear radins. On the contrary, the tails of the distribntions obtained with onr 
method have the proper exponential behavior, since they are prodnced by the solntion 
of the Schrodinger eqnation. We have performed hts with model independent techniqnes 
forcing the charge distribntions to have an exponential tail after a large valne of the radins. 
The qnality of these hts obtained by imposing this physical restriction is analogons to that 
obtained with onr mean held approach. 

Finally we want to remark that the mean helds we are generating have only a prag¬ 
matical origin and aim. Any attempt to connect them to microscopical theories of the 
nnclens is meaningless. 
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V; (Mel/) 

A 

R (fm) 

a (fm) 

208 

60.0 

0.3 

7.5 

0.6 


Table 1. Parameters of the Woods-Saxon potential (eq. 2.8) producing the density used 
to generate the pseudodata. 
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Nucleus 


SMH 

SMHWS 

exp. 

12 c 

lp3/2 

-15.53 

-16.00 

-15.96 


lpl/2 


-3.5 

-1.94 

WQ 

lp3/2 

-24.17 

-18.45 

-18.44 


lpl/2 

-16.7 

-11.97 

-12.11 


ld5/2 

-12.02 

-15.75 

-14.73 


2 sl/2 

-8.67 

-10.71 

-10.33 


ld3/2 

-9.07 

-8.64 

-8.30 

208 

lg7/2 

-11.88 

-12.78 

-11.51 


2dbl2 

-9.96 

-10.40 

-9.71 


lhll/2 

-9.10 

-9.22 

-9.37 


2d3/2 

-8.36 

-8.88 

-8.39 


3sl/2 

-7.65 

-8.12 

-8.03 


l/i9/2 

-4.21 

-5.21 

-3.80 


2/7/2 

-3.37 

-3.81 

-2.90 


H13/2 

-2.45 

-2.50 

-2.19 


Table 2. Comparison with the experimental single particle energies, in MeV. The ex¬ 
perimental valnes have been taken from ref. [16]. The SMH valnes have been obtained 
with the Hermite expansion of the potential presented in sect. 2. The SMHWS have been 
obtained expanding the central part of the potentiel on the Hermite fnnction basis, and 
nsing a Woods-Saxon shaped spin-orbit term. Note that in the SMH calcnlation the lpl/2 
level of ^^(7 was nnbonnd. 
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Nucleus 

SMH 

SMHWS 

FB 

Hermite 

12 c 

1.42 

1.25 

1.73 

1.11 

WQ 

4.44 

4.01 

3.55 

3.38 


1.35 

1.31 

1.38 

1.00 

208 

1.59 

1.52 

1.44 

1.39 


Table 3. Reduced for the fit to the experimental elastic cross sections obtained within 
our mean field model (SMH and SMHWS) and with the model independent expansion 
techinques using a Fourier-Bessel (FB) and a Hermite basis. 
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Figure captions 


Fig.l Results of the fit to the cross section pseudodata. The upper panel shows the charge 
distributions, the middle panel the central potentials and the lower panel the spin-orbit 
potentials. The full lines show the quantities used to generate the pseudodata from 
the Woods-Saxon potential. The dashed lines are generated htting the cross section 
with the potential constructed with eqs. (2.4) and (2.5), and the dashed dotted lines 
by the model where the spin-orbit potential is produced by the symmetric function 
around the hrst maximum of the derivative of V (r) starting from the external part of 
the nucleus (see text). 

Fig.2 Charge distributions and central potential produced htting the pseudodata with var¬ 
ious expansion terms of the potential, eq. (2.5). 

Fig.3 Charge distribution obtained htting the data of ref. [11]. The shaded areas 

contains all the charge distributions and potentials obtained by a set of 100 hts of the 
data normally sampled within the error band as discussed in ref. [4]. The panels A 
and B shows the distribution bands obtained with a Fourier Bessel expansion of the 
density. In A, 9 expansion coefficients have been used, and in B 13. In the panel C 
and D the darker areas show the density bands obtained with our method using 8 and 
14 coefficients respectively. The lighter areas, related to the right vertical scale, show 
the dispersion bands of the potential. We add also the information about the reduced 
averaged each set of hts. 

Fig.4 Results of the ht to the experimental data of ref. [10,11,13,14,15]. The uncertainty 
bands are produced by the Fourier-Bessel ht to the data, as discussed in ref. [4]. The 
full lines, usually lying within the uncertainty band, are generated by our mean held 
model. The dashed lines have been produced substituting the model spin-orbit term 
with a Woods-Saxon one (see text). 

Fig.5 Central (upper panel) and spin-orbit (lower panel) potentials produced by the ht to 
the experimental data. The full lines refer to ^^C, the dotted ones to the dashed 
ones to ^°Ca and the dashed dotted ones to ^°®Pb. 

Fig.6 Comparison between the central potentials of the previous hgure(full lines) and those 
generated htting the cross section data using a spin-orbit term of Woods-Saxon type 
(dashed lines). 

Fig.7 Charge distribution for the four nuclei considered presented in logarithmic scale. The 
results of our model are shown by the full and dotted lines. SMH indicates the model 
presented by the eqs. (2.4) and (2.5) where the potential has been expanded on a 
Hermite functions basis. SMHWS indicates the same model where the spin-orbit 
term has been substituted with a Woods-Saxon term. The dashed and dashed dotted 
lines have been obtained by a ht using the model independent procedure where the 
density distributions have been expanded on a Fourier-Bessel (FB) and Hermite basis 
respectively. 
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